Expansion of the Yang-Mills Hamiltonian in spatial derivatives and 

glueball spectrum 

Hans-Peter Pavel 



Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, Dubna, Russia 

and 

Institut fur Kernphysik, Technische Universitat Darmstadt 
q ■ D-64289 Darmstadt, Germany* 

(N 

■ Januar 19, 2010 

a 

Abstract 

A strong coupling expansion of the SU(2) Yang-Mills quantum Hamiltonian is carried out in the 
form of an expansion in the number of spatial derivatives, using the symmetric gauge e^A^, = 0. 
O ■ Introducing an infinite lattice with box length a, I obtain a systematic strong coupling expansion 

D \ of the Hamiltonian in A = g~ 2 ^, with the free part being the sum of Hamiltonians of Yang-Mills 

quantum mechanics of constant fields for each box, and interaction terms of higher and higher number 
of spatial derivatives connecting different boxes. The corresponding deviation from the free glueball 
spectrum, obtained earlier for the case of the Yang-Mills quantum mechanics of spatially constant 
fields, is calculated using perturbation theory in A. As a first step, the interacting glueball vacuum 
and the energy spectrum of the interacting spin-0 glueball are obtained to order A 2 . Its relation to 
\ the renormalisation of the coupling constant in the IR is discussed, indicating the absence of infrared 

■ fixed points. 

<N 

■ 1 Introduction 

O ' 

A very promising method for non-perturbative investigations of Yang-Mills theory has turned out to be 
the Hamiltonian approach [1], in particular the possibility to use the powerful variational method. 
I shall consider here the Yang-Mills theory of SU{2) gauge fields A^(x), defined by the action 



X 



S[A] := -\ J d 4 x F« U F^ , F% := d^A* - d v A^ + ge ahc A\A% , (1) 

invariant under both Poincare and scale transformations, and the local SU{2) gauge transformations 
U[oj(x)} = exp(iuj a T a /2) 

A^(x)r a /2 = U[co(x)] ^(x)r Q /2 + -d^j IT 1 ^)] . (2) 

The transition to the corresponding quantum theory is then carried out by exploiting the time depen- 
dence of the gauge transformations to put 

A a0 (x) = , a = 1, 2, 3 , (Weyl gauge) 

and impose canonical commutation relations on the spatial fields using the Schrodinger representation 

[H ai (x),A bj (y)] = i5 ab 5ij8(-K - y) — ► n ai (x) = -E ai (x) = -iS/SA ai (x) . 



'email: hans-peter.pavel@physik.tu-darmstadt.de 



1 



The physical states have to satisfy the system of equations 



(H — E)*5> = (Schrodinger equation) , 
G a (x)^ = (Gauss law constraints) , 



with the Hamiltonian ( B ai (A) = e ijk (djA ak + ^ge abc A bj A ck ^ ) 



H 



5A ai (x) 



+ B 2 ai (A(x)) 



and the Gauss law operators 



G a (x) = -i (5 ac di + ge abc A bi (x)) 



5A ci (x) 



, a = 1,2,3 , 



(3) 



(4) 



(5) 



which are the generators of the residual time-independent gauge transformations, commute with the 
Hamiltonian and satisfy angular momentum commutation relations 



[G (x), #]=0, [G (x), G b (y)] = ige abc G c (x)5(x - y) 
The matrix elements are 



<$i|0|$2> = J \{dA lk <D*0<D 2 . 

ik 



(6) 



2 Physical SU(2) Quantum Hamiltonian in the symmetric gauge 

In order to calculate the eigenstates and their energies, it is useful to implement the non-Abelian Gauss 
law constraints into the Schrodinger equation by further fixing the gauge using the remaining time- 
independent gauge transformations. One possibility, well suited for the high energy sector of the theory, 
is to impose the Coulomb gauge Xa(A) = d{A a i = describing the dynamics in terms of physical colored 
transverse gluons. I shall here choose the symmetric gauge [2] 1 

Xi{A) = e ijk A jk = (symmetric gauge) . (7) 

In contrast to the Coulomb gauge, the symmetric gauge allows for an expansion of the physical Hamilto- 
nian in spatial derivatives, which makes it very suited for the study of the infrared sector of Yang-Mills 
theory. The physical degrees of freedom in the symmetric gauge are the six components of a colorless 
local symmetric tensor field. 

The symmetric gauge corresponds to the point transformation to the new set of adapted coordinates, 
the three qj (j = 1, 2, 3) and the six elements S{ k = S k i (i, k = 1, 2, 3) of the positive definite 2 symmetric 
3x3 matrix S 

A al (q, S) = O ak (q) S kl - ^e abc ((D (q) t O T (q))^ , (8) 

where 0(q) is an orthogonal 3x3 matrix parametrized by the qi. After the above coordinate transfor- 
mation (8), the non-Abelian Gauss law constraints become the Abelian conditions 

G a $ = O — $ = (Abelianisation), 
Sqi 

that the physical states should depend only on the physical variables Si k , and the system (3) reduces to 
the unconstrained Schrodinger equation 

H(S,P)$(S) = E$(S) . (9) 

It has been proven in [3], that the symmetric gauge exists (at least for strong coupling), by showing that any time- 
independent gauge field can be carried over uniquely into the symmetric gauge. 

2 In the infinite coupling limit this transformation reduces to the polar decomposition, in which the symmetric matrix 
can be chosen to be positive definite. 
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The correctly ordered physical quantum Hamiltonian [1] in the symmetric gauge in terms of the 
physical variables -S'ifc(x) and the corresponding canonically conjugate momenta -Pjfc(x) = —iS/SSikfa) 
reads 

H(S,P) = d 3 x J d 3 y P mn (x)J/C mn | st (x,y)P st (y) + ±J d 3 X (B at (S)) 2 , (10) 

with the kernel 

/C mnN (x,y) := 6 ms 6 nt 6(x - y) - 2(x n\D m (S) *D- 2 (S) D s (S)\y t) , (11) 

the Jacobian 

J = det\*D\ , (12) 

the covariant derivative 

Di{S) k i = 5 k idi - ge k i m S m i , 

the Faddeev-Popov (FP) operator 

*D M (S) = e k miDi{S) m i = e kH di - g-fki(S) , Jm(S) = S M - 5 k itrS , (13) 
and the Green function 

(x a\*D- 2 (S)\y b) =*Drf{S)™[6{x - y)] • 
The matrix element of a physical operator O is given by 

(*'|0|*) oc f H[dS(x.)]jV'*[S]CM[S} . (14) 

A great advantage of the symmetric gauge - in contrast for example to the Coulomb gauge-, is that the 
corresponding FP operator, and hence the non-local terms of the physical Hamiltonian, can be expanded 
in the number of spatial derivatives. The Green function 

(x k\*D-\S)\y I) =*D^(S) (x) [*(x-y)] , 
corresponding to the FP operator (13), can be expanded in the number of spatial derivatives 

(x k\*D-\S)\y I) = -l 7fe Y(x)%-y) + l 7fe -\x) eabc 5W[ 7b 7 1 (x) ( 5(x-y) 



-l 7fc - 1 (x) ea6c aW[ 76 - 1 (x)^ s 5W[ 7j ] 1 (x) ( 5(x - y) 



+ ... . 



3 Expansion of the Hamiltonian in spatial derivatives 

In order to perform a consistent expansion of the physical Hamiltonian in spatial derivatives, also the 
non- locality in the Jacobian J has to be taken into account. This will be achieved in the following way. 
Writing the FP operator in the form 

*D M (S) == -gjkm(S) 8 ml - -j~l l (S)e n i i d i = -gjkm(S)*D ml (S) , 

the Jacobian J factorizes 

J = JoJ, (15) 

with the local 

Jo = det | 7 | = Y[det | 7 (x)| , det | 7 (x)| = J] (<M X ) + <M X )) . = eigenvalues of S) (16) 

x i<j 
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and the non-local J = det \*D\ . Now I include the non-local part of the measure into the wave functional 

$(S) := J~ l 'H{S) , 

leading to the corresponding transformed Hamiltonian H := J^^HJ^ 1 ! 2 , being Hermitean with respect 
to the local measure Jq 

H(S,P) = ^Q-y d 3 x J d 3 y P mn (x) j7o/C m „ N (x,y)P st (y) + lj d 3 x(B m (S)) 2 + V mcas (S) , (17) 
on the cost of extra terms 3 V meSLS from the non-local factor J of the original measure J 
V meas {S) = \Jq X J^ ss (x ) Jojd 3 ylC mnlst (x,y)A st (y) +^jd 3 xjd 3 yA mn (x.)lC mnlst (x, y)A st (y),(18) 
with the original kernel K in (11) and 

A mw (x):= ^ ln ^ = -9 (<x mrZJ-V n) - <W* fc^'V fc)) - ^(0)( 7 ~^(x) - <5 mn tr 7 - 1 (x)) . (19) 
The matrix element (14) of a physical operator O becomes the product of local matrix elements 

<*'|0|*> oc /n[^(x)n^(x) + ^(x))]$'*[S]0*[S] . (20) 

x i<j 

The transformed physical Hamiltonian (17) can be expanded in the number of spatial derivatives 

H = H + £ + ^) + £ + ... , (21) 

a \ p 7 J 

with the free part Hq containing no spatial derivatives, the interaction parts containing one spatial 
derivative, and vjf^ , Vy dd ^ A ^ containing two spatial derivatives, and so on. 

3.1 The free part H 

The free part Hq containing no spatial derivatives reads 



H = Jd 3 K- {P mn f - iS(0)[y-\(S) - 5 mn tr( 7 - 1 (5))]p mn + \^ n St^r + y (tr 2 S 2 - trS 4 ) 



,(22) 



with the spin densities 5j Spm = 2eij k Sj a P ak , z = 1,2,3 (note the factor 2). 

In order to achieve a more transparent form for the reduced Yang-Mills system I shall limit myself in 
this work to the principle orbit configurations 

< fa < fa < fa < oo , (23) 

for the eigenvalues fa , fa , fa > of the positive definite symmetric matrix S (not considering singular 
orbits where two or more eigenvalues coincide) and perform a principal-axes transformation 

S = R(a, (3, 7 ) diag (fa, fa, fa) R T (a, (3, 7 ) , (24) 

with the SO(3) matrix R parametrized by the three Euler angles x = ( a >/3>7). The Jacobian of (24) is 
\dS/d(a, (3, 7, fa\ oc sin /?]!;<.,• (4>i ~ 4>j)- The original physical variables can then be written in terms of 
the new canonical variables as (using Clebsch-Gordan coefficients) 4 

S lk = C!f lk (faff + -Ls ik (fap , P ik =Cl? lk {^f2+ ^S ik (7r 3 f ) (25) 



3 Although in principle, Kncas is part of the electric term of the Hamiltonian, I shall treat it separately in this work as 
" measure term" . 

4 For spin-1 fields sf > I use the Cartesian combinations = SfJ := (S$ - S^\)/y/2 = ef , Sf > = := i(S$ + 
S<lI)/V2 = e%S^\ = = ef , such that e.g. C\t xk := e?e?Cl£ ip , For spin-2 fields I use correspondingly the 
real combinations 4+ := (S i +l+S < L ) 2 )/V2, S?2 := -i(S$ -S ( *l)/V2, := i(Sfl +S { *l)/y/2, S[ 2 } := -(S™ -S ( ^)/V2. 
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with the spin-0 and spin-2 fields (using Wigner D-functions) 

(0 3 ) (O) := (0i + 2 + 3 ) /V3 



(03)' 



( 2 ) - 



^3 - \ (01 + 02)) Df Q { X ) + ^ (0! - 2 ) dS+(x) 



and (using 7Tj = —iS/Sfy) 

(vr 3 ) (0) := (vri + vr 2 + vr 3 ) /y/3 



(^3,6) ( 1 ) 



+ 



VT3 " i (TTi + 7T 2 )) ^ ( X ) + ^ fa - 7T 2 ) ^ + ( X ) 

n( 2 ) < \ & i n( 2 ) / \ ^ 2 i n (2) / \ ^3 
D A1+(X)— + D AlM)— —7 - + ^A2-(X)-— 



(26) 

(27) 
(28) 

(29) 



^2 - 93 93 ~ 91 

with the intrinsic spin angular momentum densities £j(x) = — i?^(x(x))5j pm (x), 

[5f in (x),e,(y)] = , [&(x),&(y)] = -ie tjk 5 3 (* - y)&(x) . 
The spin vectors 5^ pin , finally, can be written as 

Hence, in terms of the principal-axes variables, the part Hq of the physical Hamiltonian, containing no 
spatial derivatives, reads 



cyclic 



2 2id(0) , , ^ r, 

VT,- r75 -(0j7Tj - 0fc7Tfc) + 4, 



>j of"' 



The matrix elements of a physical operator O are given as 
(¥'|0|¥) oc]J y*da(x) sin /3d/3(x)d 7 (x 



2 + 0i 



,2-2 -2 



'(02 ,,2)2 ^> 



(30) 



cyclic 

[ J] d0,(x)(0 2 ( X )-02(x)) 

O<0i(x)<<fe(x)<<}!>3(x) t,j ' k 



^'*0^ . (31) 



3.2 First and second order interaction terms 



The interaction parts of first and second order in the number of spatial derivatives, needed in this work, 
can be written in the general form (A = d% + dy + c* 2 ) 



In particular, the first order magnetic part reads 

vi% n = gyJ\c\ k A2B J d 3 x (0!0 2 g id k (hf^ , 

and the second order magnetic part is 



ei=-i/A[(*fA(*f+(*SA( 



y(0A0) , y(2A2) 
y magn ' y magn 



(32) 
(33) 

(34) 
(35) 



From the second order term (35), we shall need in this work only the part Vmagn^ and the expression 



y-(S) y (S) 
M M 



1 

magn 3 



(0 3 f (0 3 f + (0 3 fj (0 3 fj 



0) 



(36) 



The first order electric term consists of transitions from spin-0 and spin-2 to spin-1 fields and therefore 
does not contribute. Of the second order electric term I shall here only need (0, = 4>j + 4>k , k cyclic) 



V 



(0A0) 



elec 



with the functions (i,j,k cyclic 



-(r-]"M'.r--(»(i-a]"(i 



Mi 



(37) 



Ti = -tt 2 + i5(0)4>i 



TTi = TTi + i8{Q)(j)i 



■ + 



% <t>l~<f> 
+ 0FWy 



7Ti + 



j + ti ,2 , 1 



and 



y(S) X (S) 
M M 



V0§ 7 4 



1 f 1 f 1 ,2 , i 7 



0f 



<5(0) 2 



\(o) 



21 



1 / 1 



f 0) 1 1 

' 2 010203 \03/ 1 2 3 

,(°\i 



^0| 



+ 11 ^ - + 



l 



\(0) 

+ 3- 



1 



5l<P2<P3 



The first order measure part, containing spin-2 and spin-3 fields 5 , reads 

5(0) [5 „ 1t /•., r[ / i \ (2) / s 3 /i i \\ (2) / 1 / 



"meas 



2A 2B 



d 3 X 



+ 



+ 3 



8g 



3M 25 



/d 3 X 



H<P2/ A > 

£ + §1 + 

*2 2 A 

/ , \ (3) 



i 



\(2) 



v 0102 \01 02 



'/A 

k 0i /A 



i 



(38) 



\(2) 



id k — 



1 



,(2) 



?i 2 JA 



+ 



M 



7s + 01. 



0l 



(3) 



M 



i<9 fc 



1 



'3/B 
\(2) 



i5i 



,0102, 



9102/ B 
(39) 



Of the second order measure term, which is very complicated, I shall need here only 



meas 



1 

+ 5 



Pi J 
\(o) 



^_ (7 + 7 

>102 \01 02 

77 Ju 



V0102 



-4 



1 



\(0) 

+ 



102. 
1 



)192 



A 



\(o) 



1 



,010203 . 



\(0) 



A - +7 



1 



\(o) 



2 s - 



\(0) 



03 



\(°h 



0102 V01 02 7 \0102 1 01 02 



and 



+ 



M M 



h4>3, 



\(0) 

A 



v 0102 \01 02 



' + if77i7' 



A 

,(o) 



x(0) 

)- 2 

/ 03 



1 



»102/ 
\(0) 



.010203 
\(0) / , 



\(0) 



,0102, 



<5(0) 5 



d 3 x 



mcas 



1 



)7a(l(i + f 

' V03 V01 02, 



\(°h 



8^3 <7 2 

5 Using the notation (</> 3 • V"3)m : = [(02^3 - <i>zi>i) + {(j>3ipi - <pii>z) + (4>iip2 ~ 4>ii\>\)\ ^ 



M 2\X-) M -2 



(40) 

(41) 

(X)) /(2>/3) 
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4 Coarse graining and strong coupling expansion in A = g 2 / 3 

I now set an ultraviolet cutoff a by introducing an infinite spatial lattice of granulas G(n, a), here cubes 
of length a, situated at sites x = an (n = (m, n,2, n-3) G Z 3 ), and considering the averaged variables 



0(n) := — f dx 0(3 

a- 3 JG(n,a) 



(where in particular <5(0) — >■ 1/a 3 ), and the discretized first and second spatial derivatives (s=l,2,3), 



A' 



1 



d s 4>(n) := lim V wjv(n)- — (0(n + ne s ) - (p(n - ne s )) 

A"— >oo — 



n=l 



2na 



<9s<K n ) : = Jim Vu;jv(n)— — ^(0(n + ne s ) + 0(n-ne s )-20(n)) 

n=l v 7 



(42) 

(43) 
(44) 

(45) 



with the unit lattice vectors ei = (1, 0, 0), &2 = (0, 1, 0), e3 = (0, 0, 1) and the distribution 

(-l)«+i(AT!) 2 

:=2 (jV- n)!(Ar + n)! ' 1 " n " " ' X>"(n) = 1 . 

The values of <9 s 0(n) and 9 2 </>(n) in (43) and (44) for a given site n and direction, say s = 1, are 
chosen to coincide with the first and second derivative, I' 2N {ani)\ n2>n3 and I'2 N {ani)\ n2>n3 6 ,respectively, 
of the interpolation polynomial l2N(xi)\n 2 ,n 3 in the x\ coordinate, which is uniquely determined by 
the series of values 4>(n\ + 71,122,113) (n = —N, ..,N) obtained via the averaging (42), and then taking 
the limit TV — > 00. Note, that the (N = 1) choice, d s (f>(n)\N=i = (<Xn + e s ) — <fi(n — e s )) /(2a) and 
dg(j)(n)\]y = i = (0(n + e s ) + 0(n — e s ) — 2</>(n)) /a 2 , which includes only the nearest neighbors n ± e s , 
would lead to the same results as (43) and (44) for the soft components of the original field </>(x), which 
vary only slightly over several lattice sites, but lead to values falling off faster than (43) and (44) for 
higher momentum components approaching n/a. 

Applying furthermore the rescaling transformation (afterwards again dropping the primes) 



4>i 



7T,; 



a a 
I obtain the expansion of the Hamiltonian in A = g~ 2 ^ 

,2/3 



2 n i ' 



H 



9 



n + XY: V<°> + A 2 I E Vf ] + £ Vf* A ) I + 0(A 3 ) 



(46) 



(47) 



with the "free" Hamiltonian 

cyclic 

^0 = E " 



2i 



(^■(n)TTj(n) - 0fc(n)7r fc (n)) 



ijk 



+C 2 (n) 



ffin)+0g(n) 
2 (n)-0 2 (n)y 



+ 2 (n)0 2 (n) 



(48) 



which is the sum of the Hamiltonians of SU (2)-Yang-Mills quantum mechanics of constant fields in each 
box, and the interaction parts, relating different boxes, 

N 



V< a >= lim 5>jv( 

W— >oo ' 
n=l 



2^ 



A? 



n=l 



C% Ml S 2 M 2 E & ^ n )] {yful [<K* + ne.)] - 3><g [0(n - ne.)]) 

n 

4E45W n )l(4W n + «e s )]+^g[0(n - ne s )]- ^^(n)]) 



(49) 



(50) 



differentiating the Lagrange interpolation polynomials I 2 n(x) with given values y n at the equidistant points x n = xo+na, 
(n = —N, —N + 1, .., N — 1, AT), once/twice at the central point xo, one obtains: I 2N (xo) = 2~D n= i w N{n)(y n — j/-„)/(2na) 
and I 2N (xq) — ^^ =1 wjv(n)(j/ n + t/_„ — 2y a )/{na) 2 with the distribution (45). For N=l, in particular, one has ^2(^0) = 
(yi -J/-i)/(2o) and 7 2 '(a;o) = (yi + 3/-1 -2y )/a 2 
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with the dimensionless and coupling constant independent terms X,y, obtained from the X, Y in (32) 
and (33) by putting X[4>] := X[0]| a=1)g=1)5(o)=1 and y[4>] := Y[</>]| a=liS=M(0)=1 . 

The expansion of the Hamiltonian in terms of the number of spatial derivatives is therefore equivalent 
to a strong coupling expansion in A = <?~ 2 / 3 . It is the analogon of the weak coupling expansion in 



^2/3 £ or sma u boxes by Liischer and Miinster[4],[5] 7 , and supplies a useful alternative to strong coupling 
expansions based on the Wilson-loop gauge invariant variables, which had been carried out by Kogut, 
Sinclair, and Susskind [6] for a 3-dimensional spatial lattice in the Hamiltonian formalism, yielding an 
expansion in 1/g 4 , and by Minister [7] for a 4-dimensional space-time lattice. 

The low energy spectrum and eigenstates of H® M at each site n appearing in (48), 



< M (n)|<On 



' {s) ( n \\$ {s) \ 



(51) 



characterised by the quantum numbers of spin S, M, have been obtained in [5], [8], [9] with high accuracy. 
It is important to note (see [9] for details), that at strong coupling, due to the positivity of the range 
(23), all states should satisfy either the (+) b.c. d ( j > p(4>)\ ( j >1= o = 0, or the (— ) b.c. <&(</>) |<£ 1= o = 
accordance with (23) and the invariance of the Hamiltonian "Ho under parity transformation <p — > - 
The spectrum is purely discrete in both cases and the lowest energies are 



m 



4.1167 , en 



The energies (relative to eo) 



(S)+ 



AS)+ 



(S)- 



8.7867 



^0 > 



(52) 



(53) 



of the lowest states for spin-0,2,3 and 4 for (+) and (— ) b.c. are summarized in Table la and lb. Spin-1 
states are absent for both cases. The underlined values correspond to stable excitations below threshold 



(= 2^ 2)+ 



), 



A*Th = 5-089 (= 2/4 2) ) , 



Mth = 3-796 

(2) 

for decay into two spin-2 excitations /j\ } (lightest in the spectrum) 



(54) 





S = 


S = 2 


S = 3 


S = 4 




(S)- 


S = 


S = 2 


S = 3 


S = 4 


i = 1 


2.270 


1.898 


8.009 


3.61 




i = 1 


3.268 


2.545 


9.250 


4.93 


i = 2 


3.857 


3.704 


10.815 


5.23 




i = 2 


5.233 


5.212 


12.78 


7.37 


i = 3 


5.09 


5.22 


13.1 


6.9 




i = 3 


6.803 


6.612 


15.38 


9.6 



(S) 

Table la and lb: Results for the first three excitation energies fi\ for(+) and (-) b.c. The under- 
lined values correspond to stable excitations below threshold (54) . The numerical errors (estimated from 
the deviation from the virial theorem, see [9]) are smaller than the last digit in the numbers given. 



5 Perturbation theory in A = g 2 / 3 

5.1 Free many-glueball states 

The eigenstates of the free Hamiltonian 



2/3 

Ho = 9 —^ M (n) 



are free many-glueball states (completely decoupled granulas). The free glueball vacuum is 

.2/3 



|0) = (g)|$o)n -> E^=Me 



9 



7 Integrating out all higher modes in a small box of size a, a weak coupling expansion for energies of the constant fields, 
E = i y]^°- efeA fc , A = [g(AMsa)] 2 ^ 3 is obtained, with the standard running coupling constant in the MS scheme. 
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(J\f total number of granulas) with all granulas in the lowest state of energy eo- The free one-glueball 
states, which in this work I choose to be momentum eigenstates, are 



\S,M,i,k) = J2e ia]l - n 

n 

the free two-glueball states, 



I* 



S>n (g) |*0>m 
m^n 



E^ ec (k)=^>»— + E; 



(S)9 



2/3 



nfree 



\(S 1 ,M 1 ,i 1 ,n 1 ),(S 2 ,M 2 ,i 2 ,n 2 )) = ® |<l>n 2 



(g) l $ 0)m 
m^ni,n 2 



p (Si,S 2 )free , (Si) 



(Si) , „(S 2 )x£? 



2/3 



/V > — • /; - 



free 
vac ' 



and so on. Matrix elements between these free glueball states are calculated using the measure (31). 



5.2 Interacting glueball vacuum 

The energy of the interacting glueball vacuum up to A 2 



-EVac — N 



2/3 



£0 + ^ B 0|vno)-A^ E E 

/3 \a,a' 11,12 Mi x + A l i 2 



+ 



(0|V^ 3 |2 n 3 t2 )(2 n 3 t2 |V^|0) ^ 

U,12 P n ^ P l2 / 



,2/3 



e + c A 2 + 0(A 3 



(55) 



is obtained using first and second order perturbation theory. 



For any of (50) I obtain, using limjv^co J2n=i( w N( n )/ n ' 2 ) = C(2) = 7I " 2 /6 , 



A? 



CO 



1st ord 



= (0|v( A) |0) = vr 2 



<*o| (4%4%) - (*o|^ 0) |*o>(*o|A:j 0) |*o> 



(56) 



For example, for the magnetic potential Vnmgn, corresponding to the Vmagn in (35), Eq. (56) becomes 



(A) 



CQ 



1st ord 7r 



magn 



(<f> m + 4 + 4>i\^o) - g(*o|& + 02 + <fe|$o> 2 J = 4.560 (3.514) . 



Here (and in the following paragraphs) the number without brackets corresponds to the (+) b.c. and 
the number in brackets to the (— ) b.c. The numerical errors are smaller than the last digit in the 
numbers given. Together with the corresponding contributions 42.323 (13.229) for the electric and 
— 16.408 (—1.782) for the measure terms V^L, I find the total first order 



co 



1st ord 



tot 



30.474 (14.962) . 



Using vf^ in (49) and limjv^oo J2n=i(. w N( n ) f 71 " 2 ) = ^C(2) = 27r 2 /3, I obtain for the contribution due 
to the vacuum polarization into a virtual pair of spin-2 particles , 



co 



2nd ord 



tt 2 v ^o||^ 2) ||< 2) )^o||yi 2) ||< 2) ) 

25 ^ „( 2 ) . „( 2 ) 

2-2, a, a' U,«2 Pii > ri 2 



+<*2 ) iiyS ) n*o><*i?iiyS ) ii*o> 



(57) 



and similarly that due to the vacuum polarization into a spin-2 and a spin-3 particle, 
2ndord ■> ^$ ||3;(3)||^( 3 )\/^„in;(2)iifls(2)\ 
35 ' 



co 



2-3 



U,»2 



Mil ' H12 
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The leading contribution to (57) comes from the Vmign — Vmagn vacuum polarization (see V^agn in (34)) 



,(9) 



r(d) 



CO 



2nd ord 



7T 2 <<&q|| (010 2 ) (2) ||< 2) X$o|| (^ 3 ) (2) 



magn— magn ii,«2 



Mil ^ Mj 2 



K 2) ||(^i^) (2) ||*o)K 2) ||(^) (2) H*o) 



+ (^ (2) ||(01^) (2) ||<J>O>(^ {2) ||(03) (2) ||^O) 



-0.399 (-0.341) 



Together with the smaller contributions -0.1516 (-0.0186) from Vnfign - v£L, -0.0295 (-0.0004) from 



Vnfias - VnSas, and the negligibly small -3.5 x 10~ 5 (-8.6 x 10~ 5 ) from V^ 3 - V£L ; 3 , I find the total 
second order 



,(9) 



,(d) ,,(d) 



2nd ord 



tot 



-0.580 (-0.360) . 



Hence 1st and 2nd order perturbation theory together give the result 



,.2/3 



4.1167 + 29.894A 2 + C(A 3 ) 



7 2/3 



8.7867 + 14.602A 2 + 0(A 3 ) 



(58) 



for the energy of the interacting glueball vacuum up to A 2 , for the (+) and (— ) boundary conditions, 
respectively. The results are summarized in Table 2. 



vacuum 


eo 


c o 


J2nd) 
c 


L C0 


co/eo 


(+) 


4.1167 


30.474 


-0.580 


29.894 


7.262 


(") 


8.7867 


14.962 


-0.360 


14.602 


1.662 



Table 2: Results for the interacting glueball vacuum for (+) and (-) b.c. The numerical errors are 
smaller than the last digits in the numbers shown. 



5.3 Interacting Spin-0 glueballs 

Including interactions V A and V 9 using 1st and 2nd order perturbation theory, we obtain the following 
energy of the interacting spin-0 glueball up to A 2 , 



7 2/3 



^+A 2 E(0^|V^|0zfc)-A 2 ( EE 



,(A) 



(0^|vg_ 2 |2 n 2 l2 )(2 n 2 J2 |V^ 2 _ 2 |0^) 



(d) 



• a,a' 11,12 
)( d ) la. q. \/a. q. 



(2) , (2) (0) 



x ^ (0ifc|Vr 3 l^3, 2 )(2 n 3 t2 |V 2 ^|0ifc) ^ 

+ m m wo I + U \ A ) 



n,i2 



J 2 ) + „(3)_ M (0) 

Ml! * H-l 2 H'l 



,2/3 



mS 0) + A 2 (cf } + cf Vfc 2 + 0((a 2 /t 2 ) 2 )) + 0(A 3 ) 



(59) 



All spin-0 glueball excitations are unstable at tree-level, except for the lowest /xf \ which is below thresh- 
old (54) for decay into two spin-2 glueballs . 

,(A) 



,(0) 



For a potential term of the general form Vg of (50), I find in first order perturbation theory 

1st ord 



M 0) i(^>S) i*i 0) > - <*oi i*-) 

- K^f l*f } > - <*o|*j 0) |*o>) <*o|*j 0) |*o> - <*o|*j 0) |*o> - (*o\4 0) \*o)) 

(60) 



-(^\x^\m^o\4 W) - <* |^ 0) |*i 0) >(*i 0) |A:f|*o> 



.(0) 



1st ord 



(^\X^\^ \4°W) + ($ |^ 0) |$i 0) )($i 0) |^ 0) |$0) 



(61) 
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For example, for the magnetic potential Vmagn corresponding to the Knagn in (35), Equs. (60) and (61) 
give 

T 2 



c (0) 



1st ord 7T 



magn 



+ 4 + ^3l^ 0) ) - <$0|# + 02 + 03l^O> ) " =K*O|01 + 02 + 3 |*i' 



2 i J.2 i /2i 



s(0)\|2 



-2/3 ((^fVl + 02 + 3 |$ 1 ° ) >($ O |01 + 02 + 3 |$O> " <$o|01 + 02 + 3 |$O> 2 ) 

= 5.296 (2.710) , 



d°) 



1st ord 



d(*o|0i + 02 + 3 |^ O) >| 2 = 0.050 (0.048) . 

magn y 



Together with the electric contributions c\ = 14.161 (4.660), c\ = 0.3977 (0.1778) and the measure 
contributions c\ = —1.813 (— 0.3597), ci = 0.0119 (0.0016), I obtain the total value from first order 
perturbation theory 



,(o) 



1st ord 



tot 



17.643 (7.011) , c5 0) " " = 0.460 (0.228) . 



1st ord 



tot 



Furthermore, second order perturbation theory leads to a change in the mass to to its virtual decay 
into two spin-2 particles or into one spin-2 and one spin-3 particle. Using in (49) I obtain for the 
case of the virtual decay into two spin-2 particles, 



2nd ord 



25 ^ ,.(2) , „(2) (0) 
2-2,a,a' U,i2 M^ + Mj 2 - Ml 



<*0||3^£ 2 >||<&^ > ><<&i 0> ||3^>||<&g ) > <*i o:> ll^£ B ^||*5? :> ><<&o||3^£»>||*g :> > 



<*£Wii*o>K 2) ii^ 



( 2 )ii-»;( 2 )ii*(°)\ jl / ( f,( 2 )ii^( 2 )ii ( T,(°)\/^( 2 )ii^( 2 )i 



(62) 



2nd ord 



— £ 



100 ^ „( 2 ) , „( 2 ) „(°) 

2-2,a,a' ii,i 2 Mil + Mj 2 ~ Ml L 



($o||^||^ 2) )(^ 0) ||^||< 2) ) 



(2)\ 



x( <*J? > ll^2 > ll*£ [>) ><^g > ll3^ > ll*o> -H <^^H3^ J 'll^o><^r''ll3^''ll*i' u ' J '> 



( 2 )ii^( 2 )i 



>( 2 )in;( 2 )iins(o)N 



-K<<^i 0) 1 13>4 2) I l^i? > X^>o| |3^4 2 > 1 1 ( <^i? ) 1 13>^> 1 1 ^>i°^ > <<^i? > 1 | |^o> 



>( 2 ) 



( 2 )||&( 2 )\ f /*( 2 )|iv( 2 )|i*(°)\/*( 2 )|iv( 2 )| 



+(^11^11*0X^11^11*? 



, (63) 



and similar expressions for the virtual decay into one spin-2 and one spin-3 particle. The leading con- 
tribution to (62) and (63) comes from Vnfign - Vnfign with c x = -3.019 (-2.5238), ci = 0.0226 (0.0194). 
Together with the smaller contributions a = -0.9947 (-1.001), ci = 0.0055 (0.0005) from Vnfign - Vnfias, 
the contributions a = -0.1189 (-0.0013), ci = 0.0036 (3 x 10~ 6 ) from V$h> - VmeL, and the negligibly 
small contributions ci = -8.3 x 10~ 5 (-2.4 x 10~ 4 ),ci = 1.2 x 10~ 7 (2.3 x 10~ 7 ) from vf\ - vf\, I 
obtain the total value from second order perturbation theory 



,(°) 



2nd ord 



tot 



-4.133 (-2.626) , c{ 



d°) 



2nd ord 



tot 



0.028 (0.020) . 



First and second order perturbation theory give the results (up to A 2 ) 



E^-(k)-E- a( 



2.270 + 13.510A 2 + C(A 3 



3.268 + 4.385A 2 + C(A 3 
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2/3 



2,„2x2n 



+ 0.488^-*; + 0((a 2 k 
a g z l 6 

2/3 

+ 0.248^fc 2 + 0((a 2 k 2 ) 2 ) 



(64) 
(65) 
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for the energy spectrum of the interacting spin-0 glueball for the (+) and (— ) boundary conditions, re- 
spectively. The results are summarized in Table 3 8 . 



spin — 


Mi 


i lst ) 
c i 


J2nd) 
c l 


Cl 


ci/mi 


-(1st) 

c l 


~<2nd) 

c i 


Cl 


l/(2/xi) 


(+) 


2.270 


17.643 


-4.133 


13.510 


5.953 


0.460 


0.028 


0.488 


0.220 


(") 


3.268 


7.011 


-2.626 


4.385 


1.342 


0.228 


0.020 


0.248 


0.153 



Table 3: Results for the interacting spin-0 glueball for(+) and (-) b.c. The numerical errors are smaller 
than the last digits in the numbers shown. 



5.4 Discussion of the results 

First I would like to comment on the relation between the glueball mass and coupling constant renor- 
malisation in the IR. Consider the physical mass 



2/3 




i H + cg 4/3 . (66) 

Demanding its independence of box size a, one obtains 

d 3 fj, + cff^ 4/3 

7(fld) = a T 9o(a) = - 9o /3 

A* c i/o 

which vanishes for the two cases, go = or g^ 3 = —c/ fx. The first solution corresponds to the perturbative 
fixed point, and the second, if it exists (c < 0), to an infrared fixed point. My result for / ^ = 
5.95(1.34) suggests, that no infrared fixed points exist, in accordance with the corresponding result of 
Wilsonian lattice QCD 9 . Solving the above equation (66) for positive (c > 0) I obtain 

5o 2/3 (Ma) = ^ + y^j -J- , a>a c :=2^I/M (67) 

with the physical glueball mass M. The minimal lattice sizes are Ma c = 11.08 (7.57), corresponding to 
a critical coupling g^\ c = 14.52 (1.55). Taking a typical physical glueball mass of M ~ 1.6 GeV [11], I 
obtain 

for M ~ 1.6 GeV : a c ~ 1.4 fm (0.9 fm) , 

which seems reasonable. The dependence of the results on the boundary conditions imposed, (+) or (— ), 
might be seen as a prescription dependence. Of course, it will be much more effective to consider mass 
ratios, as soon as, for example, the spin-2 glueball will be calculated. It would also be interesting, to 
connect the behaviour of the glueball spectrum and the bare coupling constant (67), obtained for boxes 
of large size a, with those obtained for small boxes (see [4] , [5] and footnote 7) , in order to get information 
about the intermediate region, including the possibility of the occurrence of phase transitions. 

Furthermore, I would like to remark, that Lorentz invariance imposes the following condition on the 
coefficients cf^ in (59): 

E = Vm 2 + k 2 ~ m + — k 2 -»• cj 5) = 1/(2^ ) . 

Comparison of the last two columns of Table 3 show that my result does not satisfy this requirement 
by a factor of about 2. Of course, the glueball excitation carrying non-relativistic spin-0 considered here 



8 I would like to comment here, that, using only nearest neighbor interactions (N = 1 in (43) and (44) instead of N — > oo), 
it would lead to the same ci, but a (7r 2 /6) ~ 1.64 times smaller Ci lst ' 1 |jv=i = 10.726(4.262) and a (27r 2 /3) ~ 6.58 times smaller 
c< 2nd) jjv=i = -0.628(-0.399) and hence to a 25%(12%) smaller ci|jv=i = 10.097(3.863). Similarly for the vacuum (Table 
2), it would lead to a 38% smaller c |jv=i = 18.438(9.041). 

9 In comparison, the SU(2) result from strong coupling on the lattice [7], [10]: aM = 41og(po) + 0{g$ 2 ) — > 7(50) = 
!golog(<7( 2 ) + ... does not contain infrared fixed points. 
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as a first step, does not correspond to a relativistic particle. Hence states of total angular momentum 
J = S + L, containing spin and orbital angular momentum, similar to the quark states in the Dirac 
wave-function, should be considered, e.g. the J = state (using spherically symmetric granulation) 



|$S 0) )„ (g) |*o) 

m^n 



stable 



S,i 



M 



(S) 



0/m 



where the sum is over all excitations \x 



(S) 



For simplicity, I have considered in this work only the spin-0 excitation (a4°^ 



< /i t h j underlined in Table 1, which are stable at tree- level. 

, but of course , also 



the lowest spin-2 excitations , fi^ + 



and the lowest spin-4 excitation (fJ^ ) have to be 



(2)+ (2) 

included. Most important will certainly be the inclusion of the spin-2 state fi\ ), which is lower 

in energy than the spin-0 state considered in this work. The necessary extension of the calculation to 
spin-2 and spin-4 states and the inclusion of orbital angular momentum of the lowest excitations clearly 
goes beyond the scope of this work. 



6 Conclusions 

It has been shown in this work, how a gauge invariant formulation of Yang-Mills theory on a 3-dimensional 
spatial lattice can be obtained by replacing integrals by sums and spatial derivatives by differences. This 
has been achieved by using the symmetric gauge e^Aj^ = [2] [3], and constructing the corresponding 
physical quantum Hamiltonian of SU{2) Yang-Mills theory according to the general scheme given by 
Christ and Lee [1]. In contrast to the Coulomb gauge formulation, very suitable for the description of 
the high energy sector of the theory, the symmetric gauge quantum Hamiltonian, obtained here, is very 
suitable for the IR sector, since it can be expanded in the number of spatial derivatives. The "derivative- 
free" part of the Hamiltonian is just the sum of Hamiltonians of Yang-Mills quantum mechanics of 
constant fields for each granula (here a box of size a), with a purely discrete spectrum ("free glueballs"), 
and the terms of higher and higher number of spatial derivatives describing interactions between the 
constant fields of different granulas. This expansion has been carried out here explicitly and shown to 
be equivalent to a strong coupling expansion in A = g~ 2 ^ 3 for large box sizes a. It is the analogon to 
the weak coupling expansion in g 2 / 3 by Liischer and Minister [4] [5], applicable for small boxes. Using 
the very accurate results of Yang-Mills quantum mechanics of constant fields in a box, obtained with 
the variational method in earlier work [9], the energy spectrum of weakly interacting glueballs can be 
calculated systematically and with high accuracy, using perturbation theory in A. This offers a useful 
alternative to lattice calculations based on the Wilson- loop, including the corresponding analytic strong 
coupling expansions by Kogut, Sinclair, and Susskind [6] and Minister [7]. My result for the mass of the 
interacting spin-0 glueball up to A 2 , as a first step, confirms the result of Wilsonian lattice QCD, that no 
infrared fixed points exist. Problems are a.o. the question of Lorentz invariance of the glueball spectrum. 
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